On ultrafast magnetic flux dendrite propagation into thin superconducting films 
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We suggest a new theoretical approach describing the velocity of magnetic flux dendrite pene- 
tration into thin superconducting films. The key assumptions for this approach are based upon 
experimental observations. We treat a dendrite tip motion as a propagating flux jump instability. 
Two different regimes of dendrite propagation are found: A fast initial stage is followed by a slow 
stage, which sets in as soon as a dendrite enters into the vortex-free region. We find that the dendrite 
velocity is inversely proportional to the sample thickness. The theoretical results and experimental 
data obtained by a magneto-optic pump-probe technique are compared and excellent agreement 
between the calculations and measurements is found. 
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Magnetic flux penetration in type-II superconductors 
is successfully described by the Bean critical state model 
P| . This model assumes that the slope of the flux "hills" 
is given by /j,oj c (T, B), where the critical current den- 
sity j c (T,B) is a decreasing function of the tempera- 
ture T and field B. Bean's critical state with its spa- 
tially nonuniform flux distribution is not at equilibrium 
and under certain conditions the smooth flux penetra- 
tion process becomes unstable (see review and refer- 
ences therein). The spatial and temporal development 
of this instability depends on the sample geometry, tem- 
perature, external magnetic field, its rate of change and 
orientation, initial and boundary conditions, etc. 

Instabilities in the critical state result in flux redistri- 
bution towards the equilibrium state (spatially homoge- 
neous flux throughout the sample) and are accompanied 
by a significant heat release, which often leads to the 
superconductor-to-normal-transition. The basic instabil- 
ity observed in Bean's critical state is the flux jump in- 
stability, which was discovered already in the earlyex- 
periments on superconductors with strong pinning |2|. 

The basic physics of flux jumping can be easily illus- 
trated. Assume a perturbation of temperature or flux 
occurring in Bean's critical state. This perturbation can 
be caused by an external reason or a spontaneous fluctua- 
tion arising in the system itself. The initial perturbation 
redistributes the magnetic flux inside the superconduc- 
tor. This flux motion by itself induces an electric field 
which leads to dissipation, since the electric field does not 
only act on the Cooper pairs but also on the unpaired 
electrons. This additional dissipation results in an extra 
heating which in turn leads to an additional flux motion. 
This "loop" establishes a positive feedback driving the 
system towards the equilibrium state. The flux jump- 
ing instability exhibits itself as suddenly appearing flux 
avalanche (flux jump) and heat release l^y]. 



Spatially resolved flux front patterns of Bean's criti- 
cal state instability were first observed in Nb discs with 
thicknesses in the range of d w KT 5 m to 10~ 3 m by 
means of magneto-optic imaging 4\. Wertheimer and 
Gilchrist discovered a well defined pattern of flux den- 
drites with a width w ~ 10~ 3 m and propagation veloc- 
ity v in the interval between 5m/s and 100 m/s 0. The 
dendrites velocity depended on the disks thickness, for 
smaller d a higher v was found. 

The modern magneto-optic technique allowing to in- 
vestigate flux patterns with time resolution on the order 
of 100 ps H 0] stimulated quite a few experimental 
and theoretical studies of flux front patterns arising in 
a process of smooth flux penetration pj as well as in a 
process of critical state instability development in super- 
conducting films in a transversal magnetic field. Different 
scenarios are considered resulting in a variety of flux pat- 
terns, e.g., mag netic turbulence kinetic flux front 
roughenin g |l0j| . ma g netic micro avalanches 0, > A ux 
dendrites |l3l Il4 Ha Hil It7| , thermomagnetic fingering 
[lij l , bending of flux-antiflux interface 0, H3| , and flux 
front corrugation [jlj . 

A wealth of recent experiments convincingly demon- 
strate that a propagating dendritic flux pattern driven 
by the flux jumping instability is a general phenomenon 
typical for Bean's-type critical state jH Q El El 
124 I23 ■ Indeed, the flux dendrites were observed under 
a wide variety of conditions in superconducting films of 
Nb [H[ll E|, YBazCuaOr-a [IE HE Q , Nb 3 Sn 0, 
and MgB 2 . 

It is known, that dendrite propagation in thin films 
shows velocities up to 160km/s [15| . i.e., these velocities 
are much higher than the speed of sound. This ultrafast 
motion of flux dendrites in thin superconducting films is 
a long standing and challenging problem. 

In this letter we derive a novel equation for a dendrite 



2 





FIG. 1: Magneto-optic images of a dendritic flux pattern in a 
YBCO film with the thickness d = 330 nm subjected to a field 
of B a = 15 mT. (a) Final state (after ~ 10 s) of a dendritic flux 
pattern with superimposed current distribution shown by the 
arrows. The length of the arrows is proportional to the local 
current density, (b) The absolute value of the current density 
is shown. The bright areas indicate high current densities. 



tip velocity and demonstrate an excellent agreement be- 
tween the theoretical results and experimental data for 
the propagation velocity of a single flux dendrite branch. 

Dendritic flux structures which can be considered as 
a set of single flux branches originating from a certain 
area were observed in numerous experiments [lq . In the 
case of a dendritic structure with few branches the single 
branches do not affect each other and the propagating 
substructures can be treated as a moving flux jump in- 
stability localized at the tip of the dendrite branches. A 
typical magneto-optic image of a "dilute" dendrite in its 
final state is shown in Fig. ^a). Superimposed are the 
current streamlines as determined by an inversion scheme 
[3fj| . In Fig.QJb) the absolute value of the current density 
is shown. It is worth mentioning that the center of the 
dendrite is current free and that the current follows the 
dendrite branches. The current density decreases rapidly 
with distance from the dendritic structure. 

These experimental observations allow for a straight- 
line flux dendrite model, which we use in our calculations. 



This model assumes the following: 

(a) The current of a straight-line dendrite first flows 
parallel to the sample edge, then closely follows the con- 
tour of the dendrite branch until flowing parallel to the 
sample edge again as shown in Fig. |21 

(b) In the current carrying areas, the superconductor 
is in the flux creep regime and thus the current density 
j depends on the electric field E as 

j = j c (E/E f n , (1) 

where j c is the critical current density, n and Eq are 
the parameters characterizing the current density-electric 
field curve (at E = Eq we have j — j c ). It is common 
to define j c as the current density at Eq — 10~ 4 V/m, 
for high-T c superconductors n ~ 10 but decreases with 
the applied magnetic field [2!| . Eq. (JIJ yields the electric 
field dependent conductivity 

Next, we denote the radius of the dendrite tip as po 
and the width of the current carrying area as 

P P = B cS /p j c , (3) 

where -B c ff = -Bin — B out , B m is the field inside the den- 
drite, and -B ou t is the field outside the tip of the dendrite. 

Consider now the flux front stability at a tip of a mov- 
ing flux dendrite in the framework of the model developed 
to treat the flux jump instability near a semicircle inden- 
tation at the sample edge |26l 1271 128| . This approach is 
based on the assumption that the flux jumping instability 
develops much faster than the magnetic flux diffusion. In 
the flux creep regime of low-T c and high-T c superconduc- 
tors this assumption holds with a high accuracy |2fl. l27| . 




FIG. 2: Current lines for the straight-line magnetic flux den- 
drite. The full line to the left marks the strip edge. 
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It follows from the general approach that the stability 
margin of a flux jumping instability is determined by the 
existence of a nontrivial solution of the equation [23, 
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where 9 is the temperature perturbation, A is the heat 
conductivity, E is the electric field generated by a time 
dependent magnetic field, the parameter q is given by 



tanqd = h/Xq, 



(5) 



and h is the heat transfer coefficient to the coolant. The 
boundary condition to Eq. (@J is nV6* = at the edge 
of the film and n is the unit vector perpendicular to the 
edge of the film. It is clear from Eq. J3J that the flux 
front stability is highly sensitive to the electric field E 
generated by the varying magnetic field [2(| [2?], • 

The dendrite tip motion results in an electric field E, 
which is parallel to the current density j. We consider 
this field similar to the consideration of the electric field 
generated by a varying magnetic field at a semicircular 
indentation with a radius po m a superconducting film 
with a straight edge |2^| . This approach results in 



E w B in pl/po 



(6) 



Assuming that po < Pp we estimate the magnetic field 
rate in the vicinity of a dendrite tip as 



-Bin « vB ef f/p p 



(7) 



Combining the Eq.[S]and[7|we find that the electric field 
generated at the inner edge of a moving flux dendrite tip 
can be estimated as E ss vB c sp p / pq. 

Next, we assume a high thermal boundary resistance 
in the films, which means that hd <C A. This assumption 
can be justified using typical values A = 0.1 WK^ 1 m _1 , 
d = 330 nm, and h = (700 to 1.3 x 10 3 ) WK^m" 2 HJ 
Is^ . In this case of hd <C A we find from Eq. (JjJJ that 
the value of q 2 « h/Xd and the stab ility criterion for flux 
jumping in thin film takes the form 
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Consider now the critical current density to be linear 
in temperature, i.e., j c = j (l — T/T c ). In this case we 
use Eq. (J7J to rewrite Eq. JSJ as follows 
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where 7 ~ 1 is a numerical factor and jo is the critical 
current density at T = 0. The criterion given by Eq. © 
gives the lower limit of a dendrite tip speed. If we assume 
that po « p p and use Eq. J2J then 
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FIG. 3: The time dependence of the length of a magnetic flux 
dendrite s(i). The solid line is the solution of Eq. 112H and 
the solid dots are the experimental data [ThI |. 



We compare now the results obtained by Eq. (|10f> and 
our experimental data. To measure the time dependent 
dendrite length s = s(t) we used a magneto-optic pump- 
probe setup 0|. The dendrites where nucleated at the 
edge of a square YBCO thin film sample by focusing a 
laser beam onto the film surface. The magnetic field 
was applied prior to the laser pulse. We observed two 
qualitatively different stages of dendrite propagation. In 
the first few nanoseconds we observed an extremely high 
velocity on the order of 160km/s, later on this velocity 
decreased to a value of 18 km/s. For experimental details 
see Ref. 0] . The existence of these two distinct regions 
of dendrites propagation can be easily understood using 
Eq. (|10|l . Indeed, as long as a dendrite crosses the crit- 
ical state area the field B out is decreasing, therefore the 
value of B c ff = B m ~ B out is increasing and consequently 
the velocity of the dendrite is decreasing. After the den- 
drite tip crosses the critical state area its velocity stays 
constant as the dendrite runs in a vortex free area where 
B c g is a constant. 

The time dependence of the dendrite length s = s(t) 
can be calculated using Eq. (f 1C)|) . The effective field B c g 
is the crucial parameter for this calculation. To find B c g 
we assume that the magnetic field inside a long supercon- 
ducting strip is a good approximation for the distribution 
at the center of a square superconducting thin film |33J 
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where a and b are the half widths of the superconductor 
and of the vortex free area, and c = \Jl — b 2 /a 2 . Sub- 
stituting Eq. (|11|) into Eq. I|10|) we obtain a differential 
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FIG. 4: The dependence of the flux dendrite velocity on the 
sample thickness v(d). The dashed line is a fit of v cc 1/d 
revealing the dependence given by Eq. J7J and the solid dots 
are the experimental data [34| . 



equation for the dendrite propagation 



ds 
~dl 



a 



Bin — -B out (s) ' 



(12) 



with a = 2jhT c /ndjQ. Based on our experimental data 
we assume that the field B m is constant and by a factor 
of 1.9 larger than the applied magnetic field. 

A numerical solution of Eq. I|12fl yields the solid line 
in Fig. [3J We used for this plot the values d = 330 nm, 
a = 5 mm, b = 4.4 mm, T c = 90 K, j = 1.5 x 10 9 A/m 2 , 
h = 10 4 WK- 1 m- 2 , T c = 90K, 7 = 1, and n = 6. We 
take s(0) = a — I with t = 0.1 mm as an initial condition 
to avoid the singularity of B out (s) at t = 0. 

To check Eq. (jlUII further we compare the calculated 
velocities with the velocities obtained from line-focus 
measurements 0|- If a nne focus is applied, the den- 
drites never run in the critical state, but penetrate into 
flux free area. In this expected from Eq. (|To"|l . 

we don't find a regime with increased velocities, however 
we find a thickness dependence. In Fig. 21 one can see 
the experimentally obtained thickness dependent veloc- 
ity and a fit v oc 1 / d. One reason for the slight deviations 
between the theory and experiment may be that we had 
to use different YBCO films to obtain the data, i.e., the 
values for parameters like T c or B c h may vary. 

The main result of this letter is Eq. (|T0")l . It describes 
the dynamics of a single flux dendrite and it was shown, 
that an excellent agreement with the experimental data 
has been achieved. To the authors' knowledge this is the 



first time that a theory explains the observation of the 
fast and slow penetration velocities. 
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